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Abstract 

We present high-resolution soft-X-ray spectra of the 
prototypical Seyfert 2 galaxy, NGC 1068, taken with the 
XMM-Newton Reflection Grating Spectrometer (RGS) and 
the Chandra Low Energy Transmission Grating Spectrom- 
eter (LETGS). Its rich emission-line spectrum is domi- 
nated by recombination in a warm plasma (bright, nar- 
row radiative recombination continua provide the "smok- 
ing gun"), which is photoionized by the inferred nuclear 
power-law continuum. Radiative decay following photoex- 
citation of resonant transitions also provides an impor- 
tant contribution. A self-consistent model of a photoion- 
ized and photoexcited cone of gas is capable of repro- 
ducing the hydrogcnic/hcliumlikc ionic line series in de- 
tail. The ratio of photoexcitation to photoionization in 
the cone provides important geometric information such 
as the radial ionic column densities, which are consistent 
with absorption measurements (the warm absorber) in 
Seyfert 1 galaxies. This strongly suggests that the emis- 
sion spectrum we observe from NGC 1068 emanates from 
its warm absorber. The observed extent of the ionization- 
cone/warm-absorber in NGC 1068 of about 300 pc implies 
that a large fraction of the gas associated with generic 
warm absorbers may typically exist on the hundreds-of- 
parsec scale rather than much closer to the nucleus (e. g., 
less than a parsec). Spatially- resolved spectroscopy using 
the LETGS of two distinct emission regions yields two no- 
ticeably different spectra. We show that these differences 
are solely due to differing radial column densities. A fairly 
flat distribution in ionization parameter ^ — Lx jn^r^ 
(over at least log^ ~ 0-3) is necessary to explain the 
inferred radial ionic column densities of all spectra. We 
show that this must primarily be due to a broad den- 



sity distribution /(?^c) 



at each radius, spanning 



roughly rif. ~ 0.1-100 cni^^^. Additional contributions to 
the soft-X-ray emission from hot, coUisionally-ionizcd gas, 
if present, make a negligible contribution to the spectrum. 

Key words: galaxies: individual (NGC 1068) — galaxies: 
Seyfert — line: formation — X-rays: galaxies 



1. Introduction 

In the unified model of active galactic nuclei (AGN), the 
observational properties of a particular AGN are deter- 



mined simply by its orientation (Miller fc Antonucci 1983 



Antonucci fc Miller 1985). For Seyfert 1 galaxies, we di- 



rectly observe the intrinsic nuclear continuum, partially 
absorbed by photoionized outfiowing material (warm ab- 
sorber). We observe Seyfert 2 galaxies at an angle nearly 
perpendicular to the Seyfert 1 orientation. From this van- 
tage, the intrinsic nuclear continuum is highly obscured 
(by the "dusty torus" ) , and the spectrum is dominated by 
reprocessed emission from outlying clouds filling an ion- 
ization cone, which may be related to the warm absorber. 

The first high-resolution soft-X-ray spectrum obtained 
of a S eyfert 2 galaxy w a s the Chandra HET GS of Markar- 
ian 3 (|Sako et al. 200C| ). |Sako et al. (2000)1 interpreted the 
bulk of this spectrum as due to recombination/radiative 
cascade following photoionization. Radiative decay follow- 
ing photoexcitation is also required to explain the en- 
hanced resonance lines in the heliumlike triplets. This is 
consistent with outlying clouds irradiated by the inferred 
nuclear continuum, as predicted in the unified model Kro- 
hk fc Kriss (1995)| . 

We present below the first high-resolution X-ray spec- 
trum of the X-ray- brightest, prototypical Seyfert 2 galaxy, 
NGC 1068, which was obtained with XMM-Newton RGS 



(Kinkhabwala et al. 2002b). Using a simple model of an 
irradiated gas cone, we are able to find an excellent fit 
to the detailed soft X-ray spectrum of NGC 1068 (for all 
hydrogenic/heliumlike ions). From the similarity of our 
inferred radial column densities to directly-measured ab- 
sorption column densities (warm absorber) in Seyfert 1 
galaxies, we argue that the observed hundreds-of-parsec- 
scale ionization cone of NGC 1068 is identical to the warm 
absorber in this AGN. A su bsequent spectrum ob tained 
with the Chandra LETGS (|Brinkman et al. 200^ ) con- 
firms the RGS results, but, more importantly, allows for 
spatially-resolved spectroscopy of this extended source. 
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2. Full RGS/LETGS Spectra of NGC 1068 

Line emission dominates the RGS /LETGS soft-X-ray spec- 
tra of NGC 1068 shown in Fig. |l^ (for ease of comparison, 
all spectral plots appear at the end of the paper). Emis- 
sion lines from hydrogenic and heliumlike C, N, O, Ne, 
Mg, Si, and S are all clearly detected. Fe L-shell emis- 
sion lines from Fe XVII to Fe XXIV are present as well 
(unlabelled) , with numerous transitions scattered between 
9 A and 18 A. Many higher-order resonant transitions 
(Is^np) in the hydrogenic and heliumlike ions labelled 
(3-5 are prominent, with evidence for strong emission from 
even higher order transitions as well. Several unidentified 
features at longer wavelengths (e. g., at 27.92, 30.4, 34.0- 
34.6, and 36.38 A) are likely due to L-shell emission from 
mid-Z ions such as sulfur. The fluorescent lines of neutral 
Fe (Fe+°) and Si (Si+°) are clearly detected in the LETGS 
spectrum. We see no significant continuum emission in the 
spectrum. 

All lines appear broader than expected for a monochro- 
matic source and several lines show evidence for significant 
blue shifts. These are due to intrinsic velocity distributions 
(hundreds of km/s) and are not instrumental effects. 

The spectrum also includes very prominent radiative 
recombination continua (RRC) for hydrogenic and heli- 
umlike C, N, and O, which are produced when electrons 
recombine directly to the ground state in these highly- 
ionized ions. RRC are smeared out for hot coUisionally- 
ionized gas, but are narrow, prominent features for cooler 
photoionized gas. The narrow width of these RRC pro- 
vide a direct measure of the recombining electron temper- 



Seyfert 1 



ature (Liedahl & Paerels 1996; Liedahl 1999), which for 
NGC 1068 is fcT - 2 - 10 eV. 



3. Model of Irradiated Gas Cone 

We have constructed a fully self-consistent model of a gas 
cone (warm absorber) irradiated by a continuum source 
(Fig. p. The nuclear region comprising the black-hole/ 
accreti on-disc /comptonized-halo system is depicted as the 
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Figure 1. Schematic of our unified, warm- 
ahsorber/ ionization- cone model for Seyfert 1 and 2 
galaxies (not to scale). The labels "Seyfert 1" and 
"Seyfert 2" indicate the directions of the observer in the 
two cases. 



of fitting, we have incorporated our model into XSPEC 
( Arnaud 199(: ) as a local model ("photo"). A complete 
discussion of the atomic calculations and astrop hysical as- 



central l black circle. Obscuration by the "dusty torus" is wala et al. (2002a 



sumptions underlyi ng our model is presented in Kinkhab- 



shown (in cross section) as two clouds on either side of the 
nucleus. The parameters which comprise our two inferred 
global model parameters, the covering factor times nuclear 

lumin osity fLx and the radial velocity width ir^^'^, are much brighter than its counterpart in the SW ( Young et al 



For the specific case of NGC 1068, the intrinsic con- 
tinuum is likely to be completely obscured, contributing 
no flux to the soft X-ray regime. Also, the NE cone is 



indicated, as well as the individual radial ionic column 
densities N^^^. Lx is the total luminosity in a power-law 
with reasonable values for the index of F = —1.7 and 
energy range of 13.6 eV-100 keV. 

Photoionization and photoexcitation in the ionization 
cone produce absorption features (warm absorber) in the 
Seyfert 1 view, but in the Seyfert 2 view, the inverse pro- 
cesses of recombination/radiative cascade and radiative 
decay, respectively, produc e line emi ssion (Fig. ||). We use 
the new atomic code FAC ( |Gu 2002 ) to calculate the rele- 
vant atomic data for hydrogenic/heliumlike ions. For ease 



2001 ), therefore, the covering factor / applies to the NE 



cone alone (hence the asymmetry of Fig. 

3.1. Column Density 

We show the effect of varying the radial ionic column den- 
sity on O VII in Fig. |^. The top three panels on the left 
show the radial "Seyfert 1" view through the outflow and 
down to the nucleus for radial column densities in O VII 
of 10^^, 10^^, and 10^^ cm~^. The corresponding top three 
panels on the right show the "Seyfert 2" view roughly 
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Figure 2. Simplified Grotrian diagram of the relevant 
atomic processes in hydrogenic/heliumlike ions. Photoex- 
citations and photoionizations appear as absorption fea- 
tures in the Seyfert 1 view (left), but as line emission 
through the inverse processes of radiative decay and recom- 
bination/radiative cascade in the Seyfert 2 view (right). 



perpendicular to the axis of outflow and from which the 
nucleus is completely obscured (Fig. |]). All photons ab- 
sorbed out of the power law in the "Seyfert 1" spectrum 
are reprocessed and reemitted to generate the "Seyfert 2" 
spectrum. Radiative decay following photoexcitation dom- 
inates the "Seyfert 2" spectrum at low column densities, 
whereas recombination/radiative cascade following pho- 
toionization dominates at high column densities (for com- 
parison, pure recombination is shown in the lower right 
panel). For hydrogenic ions, the behavior is similar (omit- 
ting the intercombination/forbidden lines). 

3.2. Velocity Width 

In Fig. 1^, we show the effect of varying the radial velocity 
width for O VII. Larger a]f^'^ enhance the importance of 
photoexcitation relative to photoionization. 

4. Model Fit to the RGS Spectrum 

We present our fit to the RGS spectrum of NGC 1068 in 
Fig. |l^. The parameters used for the fit are given in the 
first column of Table |l|. To convert observed flux to lumi- 
nosity for the spectra of NGC 1068, we assume a distance 

The column 



of 14.4 Mpc ( p land-Hawthorn et al. 1997| ). 
density of neutral hydrogen in our galaxy is taken to be 
Nh ~ le20 cm^^ (using the high-quality soft-X-ray spec- 
trum itself), which is somewhat lower than the standard 
column density to this source ( Kinkhabwala et al. 2002b| ). 
For all fits in this paper, we assume recombining electron 
temperatures of kT = 2.5 eV for C V, AiT = 3 eV for 
N VI, fcT = 4 eV for C VI, N VII, O VIII, and O VIII, 
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Figure 3. Effect of differing radial column densities on the 
reemitted spectrum for O VII with the Seyfert 1 view in the 
top three panels on the left and the Seyfert 2 view in the 
corresponding top three panels on the right. The bottom 
left panel shows ionic cross section with the photoexcita- 
tion/photoionization boundary indicated. The bottom right 
panel shows the spectrum expected for pure recombination. 
All spectral axes are linear, but have arbitrary normaliza- 
tion. Throughout, we take a radial gaussian distribution 
with al'^'^ = 200 km/s, transverse velocity distribution with 
cr'^^ = 400 km/s, and temperature kT — ieV. 
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Figure 4.. The effect of different cr]f^'^ for the O VII spec- 
trum using N^^^ — 10^^ cm""^ , with the Seyfert 1 view 
on the left and the Seyfert 2 view on the right. Spectra on 
the right-hand- side were convolved with the same trans- 
verse velocity distribution a"^^ 
sume kT — Z eV. 



400 km/s; we also as- 



and kT = 10 eV for all other ions (which have poorly- 
determined RRC). 

Velocity broadening of all model emission lines by (7°'^^ = 
400 km/s is necessary to fit the observed lines. 
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Figure 5. The final fit to C VI including recombina- 
tion/radiative cascade following photoionization and ra- 
diative decay following photoexcitation (top). Recombina- 
tion alone (bottom) is unable to explain the excess emis- 
sion in all resonant lines np^ls. 
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Figure 6. The final fit to O VII including recombina- 
tion/radiative cascade following photoionization and ra- 
diative decay following photoexcitation (top). Recombina- 
tion alone (bottom) is unable to explain the excess emis- 
sion in all resonant lines np^ls. 



We illustrate the relative contributions of photoion- 
ization and photoexcitation to the ionic line series for 
C VI and O VII in Figs. || and ||, respectively. Pure re- 
combination is unable to explain the anomalous strength 
of the higher-order-series transitions. However, the self- 
consistent addition of photoexcitation allows for an excel- 
lent overall fit. An additional coUisional gas component 
instead of photoexcitation would be insufficient to explain 
the higher-order series transitions (Fig. 0). 
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Figure 7. We demonstrate how to distinguish between hot 
coUisional gas (e. g., starburst region) and photoexcita- 
tion. Starting with pure recombination following photoion- 
ization ("PI" - top two panels, kT = 4 eV), we self- 
consistently add radiative decay following photoexcitation 
("PE"- middle left panel) assuming reasonable ioniza- 
tion cone parameters, or an additional hot gas compo- 
nent in coUisional ionization equilibrium ("CIE" - mid- 
dle right panel, kT = 150 eV) to obtain the bottom two 
panels. Note that both bottom panels have similar triplet 
ratios, implying that using the triplet alone is insuffi- 
cient to discriminate between these two scenarios. How- 
ever, the "AGN ALONE" panel has significantly stronger 
higher- order- series transitions (including the RRC) than 
the "AGN-f-STARBURST" panel, demonstrating the diag- 
nostic importance of these transitions. (Normalization in 
each panel is arbitrary.) 



5. LETGS Spatially- Resolved Spectroscopy 

With the LETGS, it is possible to perform spatially-resolved 
spectroscopy by making cuts in the cross dispersion direc- 
tion. NGC 1068 provides perhaps the best example of this 
capability, since it shows evidence for two fairly separated 
emission regions, which we denote simply as "primary" 
and "secondary." The two cross-dispersion regions we use 
for the "primary" and "secondary" spots are shown in 
Figs. ^ and ^. 

Spectra of the "primary" and "secondary" spots are 
given in Figs, p^and p^ , respectively, along with their cor- 
responding best fit models. The fit parameters we obtain 
are given m Table 0. The same fLx parameter works fairly 
well for all spectra, however a slightly different radial ve- 
locity width crjj'^'^ is preferred by the "primary" spectrum. 
This illustrates the relative factor-of-two uncertainty in 
all fit parameters. However, this uncertainty does not af- 
fect any of our conclusions. All LETGS models have been 
convolved with the specific zero-order dispersion profile. 
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Table 1. Fit parameters for Figs. and [7^. 'H' and 

'He' indicate hydrogenic and heliumlike, respectively. 
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Figures. Zero-order LETGS image (logarithmic) of 
NGC 1068 oriented with N up and E to the left. There 
are two general regions which we denote as the "primary" 
(brightest spot) and "secondary" (accompanying NE spot). 
The dispersion axis is parallel to the green lines, which de- 
marcate the cross dispersion regions we used to generate 
the "primary" and "secondary" spectra. The apparent line 
of emission in the cross-dispersion direction and centered 
on the source is instrumental (due to CCD readout). The 
high-energy edge of the dispersed spectrum can be seen in 
the upper right and lower left of the image. 



The "secondary" spot, which is roughly coincident with 
the radio outflow lobes and the edge of the UV ioniza- 
tion cone, exhibits a noticeably different spectrum from 
the "primary" spot. The heliumlike resonance lines of Ne 




Figure 9. Zero-order LETGS image of NGC 1068. The rel- 
ative brightness of the "primary" and "secondary" spots 
are clearly shown. The straight lines in the "factors-of- 
two" contour image are the same as the green lines in 
Fig. 1^, denoting the extraction regions for the "primary" 
and "secondary" spectra. 



and O are relatively stronger compared to the forbidden 
lines than in the "primary" spot. This enhancement might 
be due to an additional coUisional component, perhaps 
outflow-shock-heated gas. However, this is ruled out by 
the strength of several highcr-order-series lines and the 
normalization of the hcliumlikc intercomination/forbidden 
lines to the RRC (sec Fig. for an explanation). The ob- 
served enhancements are in fact simply due to a lower 
column density through the material (Fig. ||). 



6. Conclusions 

We have shown that the XMM-Newton RGS spectrum of 
the soft-X-ray emission from NGC 1068 is due entirely 
to recombinations and radiative cascade in an ionized gas 
cone, which is photoionized and photoexcited by the in- 
ferred nuclear continuum. A simple model of a warm, 
photoionized cone is capable of explaining all hydrogenic 
and heliumlike ion spectra in detail. The values we infer 
for the radial ionic column densities are similar to col- 
umn densities observed in absorption (warm absorber) in 



Seyfert 1 galaxies ( 


Kaastra et al. 2000; 


Kaspi et al. 2000; 


Branduardi-Raymont 2001 




Sako et al. 2001 


). Since the 



ionization cone in NCC 1068 is spread over a region of 
~300 pc (Young et al. 2001), this implies that generic 



warm absorbers have typical sizes of hundreds of par- 
sec, rather than existing very close (e. g., < 1 pc) to the 
nucleus. A subsequent spectrum obtained with Chandra 
LETGS confirms the RGS results for the spectrum of the 
entire object, but also allows for spatially-resolved spec- 
troscopy of the two bright spots resolved in the zero-order 
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image. The difference in spectra between these two spots is 
due entirely to different radial ionic column densities, and 
not, for example, to an additional contribution from hot, 
outflow-shocked gas. In Fig. 10, we show that the temper- 



atures we obtain for all spectra using well-measured, non- 
blended RRC arc fairly consistent with (though slightly 
higher than) temperatures predicted from a self-consistent 
simulation using XSTAR ([Kallman fc Krolik (1995)|) of an 



optically-thin photoionized plasma. 
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Figure 11. Fractional ionic abundances for the heliumlike 
(left), hydrogenic (middle), and bare (right) charge states 
for several ions as a function of ionization parameter ^ = 
5 Lx/n^r^. 



Figure 10. Recombination emissivity peak (for formation 
of the H '- or He '-like ionic species) calculated using 
XSTAR with incident power-law spectrum (T — —1.7, 
13.6 eV< E <100 keV). The FWHM of the emissivity for 
each ion is ^ 1. The x's and o's give the confidence 
interval and best fit, respectively, for each ion temperature 
measured from the ion RRC. 



The ionic column densities we infer allow us to probe 
the ionization parameter distribution of the gas. We show 
the predicted fractional ionic abundances as a function of 
ionization parameter ^ = Lx/uer'^ in Fig. |ll|. We point 
out the relatively equal column densities inferred for the 
hydrogenic and heliumlike species of each element in each 
spectrum in Table |^. No single ionization parameter is ca- 
pable of reproducing this result. Instead, a rather flat dis- 
tribution in ionization parameter is necessary. This could 
be obtained by assuming a spatially-stratified, single-density 
ionization cone (hence varying only r) or assuming an 
intrinsic density distribution at each radius. The strik- 
ing overlap of O III and soft-X-ray emission regions in 



NGC 1068 ( [Young et al. 200lD coupled with the presence 
of a similar range of ions in the "primary" and "secondary" 
spots, which are located at different distances from the nu- 
cleus, favor the latter interpretation. The distribution in ^ 
then is mostly due to a distribution in ric of /(tt-c) cx 
over several or ders of magnitude (roughly ~ 0.1-100) 



at each radius ( Kinkhabwala et al. 2002b ; Brinkman et al 
20020 
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Figure 12. The RGS1/RGS2 (m=-l order) (blue) and LETGS (combined m=-l,+l orders) (red) spectra of NGC 1068. 
Note the overall flux and wavelength agreement between both intruments. The spectrum has been shifted to the NGC 1068 
rest frame (z — 0.00379/ The superscripts 'H' (hydrogenic) and 'He' (heliumlike) refer to the final-state ion (e. g., 
'C^°RRC' refers to recombination forming heliumlike carbon.) Unlabelled features at\<\d>A are all due to Fe L-shell 
transitions, whereas unlabelled features at \ > 18 A are due to mid-Z-element (mostly S and Si) L-shell transitions. 




Figure 13. RGSl spectrum (blue) and corresponding hydrogenic/heliumlike ion fit (red) using parameters listed in 
Table ^. The spectrum has been shifted to the NGC 1068 rest frame (z = 0.00379j. Model wavelengths are those 
expected for no excess velocity shifts. We employ the same labelling convention as in Fig. W^. 



9 




Figure 14- LETGS spectrum (m=-l/+l orders) of the "primary" spot (blue) and corresponding hydrogenic/heliumlike 
ion fit (red) using parameters listed in Table^. The spectrum has been shifted to the NGC 1068 rest frame (z — 0. 00379 j. 
Model wavelengths are those expected for no excess velocity shifts. We employ the same labelling convention as in 
Fig. |7|. 
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Figure 15. LETGS spectrum (m=-l/+l orders) of the "secondary" spot (blue) and corresponding hydrogenic/heliumlike 
ion fit (red) using parameters listed in Table^. The spectrum has been shifted to the NGG 1068 rest frame (z = 0.00379J. 
Model wavelengths are those expected for no excess velocity shifts. We employ the same labelling convention as in 
Fig. Note the enhanced resonance line in the O and Ne heliumlike triplets. Several higher-order (np-^ls, n > 2) 
lines also appear (e. g., 'O^''^' and 'N^7'j, providing clear evidence for photoexcitation rather than an additional hot, 
collisional gas component (see Fig. 



